Mixed:
data %>% mutate(met = numCs/(numTs + numCs)) -> data
data$tiles = data$poz %/% 1000
data %>% group_by(tiles, chr) %>% summarise(mean_met = mean(met), n = n()) %>%
mutate(start = tiles * 1000) -> avg
beta_mixed %>% left_join(avg, by = c("start" = "start", "chr" = "chr")) %>%
mutate(p0 = exp(mean_met)/(1+exp(mean_met)), p1 = exp(mean_met + beta)/(1+exp(mean_met+ beta))) %>%
mutate(dist = p1/p0) %>% arrange(desc(dist)) -> beta_mixed
ggplot(beta_mixed, aes(y = dist, x = mean_met)) + geom_point(size = 0.1) +
ggtitle('p1/p0 w zaleznosci od sredniej globalnej')
## Warning: Removed 564 rows containing missing values (geom_point).

ggplot(beta_mixed, aes(y = dist, x = n)) + geom_point(size = 0.1) +
ggtitle('p1/p0 w zaleznosci od pokrycia')
## Warning: Removed 564 rows containing missing values (geom_point).

ggplot(beta_mixed, aes(y = dist, x = p_value)) + geom_point(size = 0.1) +
ggtitle('p1/p0 w zaleznosci od pValue')
## Warning: Removed 564 rows containing missing values (geom_point).

ggplot(beta_mixed, aes(y = p_value, x = n)) + geom_point(size = 0.1) +
ggtitle('pValue w zaleznosci od pokrycia')
## Warning: Removed 555 rows containing missing values (geom_point).

beta_mixed %>% filter(mean_met <= 0.9 & mean_met >= 0.1) %>% arrange(desc(dist)) -> beta_mixed
Mixed z efektami losowymi
corr_mixed %>% left_join(avg, by = c("start" = "start", "chr" = "chr")) %>%
mutate(p0 = exp(mean_met)/(1+exp(mean_met)), p1 = exp(mean_met + beta)/(1+exp(mean_met+ beta))) %>%
mutate(dist = p1/p0) %>% arrange(desc(dist)) -> corr_mixed
ggplot(corr_mixed, aes(y = dist, x = mean_met)) + geom_point(size = 0.1) +
ggtitle('p1/p0 w zaleznosci od sredniej globalnej')
## Warning: Removed 3187 rows containing missing values (geom_point).

ggplot(corr_mixed, aes(y = dist, x = n)) + geom_point(size = 0.1) +
ggtitle('p1/p0 w zaleznosci od pokrycia')
## Warning: Removed 3187 rows containing missing values (geom_point).

ggplot(corr_mixed, aes(y = dist, x = p_value)) + geom_point(size = 0.1) +
ggtitle('p1/p0 w zaleznosci od pValue')
## Warning: Removed 3187 rows containing missing values (geom_point).

ggplot(corr_mixed, aes(y = p_value, x = n)) + geom_point(size = 0.1) +
ggtitle('pValue w zaleznosci od pokrycia')
## Warning: Removed 629 rows containing missing values (geom_point).

corr_mixed %>% filter(mean_met <= 0.9 & mean_met >= 0.1) %>% arrange(desc(dist)) -> corr_mixed